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ABSTRACT 



Context. Very few massive young stellar objects (MYSO) have been studied in the infrared at high angular resolution due to their 
rarity and large associated extinction. We present observations and models for one of these MYSO candidates, NGC 3603 IRS 9A. 
Aims. Our goal is to investigate with infrared interferometry the structure of IRS 9A on scales as small as 200 AU, exploiting the fact 
that a cluster of O and B stars has blown away much of the obscuring foreground dust and gas. 

Methods. Observations in the N-band were carried out with the MIDI beam combiner attached to the VLTI, providing spatial informa- 
tion on scales of about 25-95 milli-arcseconds (mas). Additional interferometric observations which probe the structure of IRS 9 A 
on larger scales were performed with an aperture mask installed in the T-ReCS instrument of Gemini South. The spectral energy 
distribution (SED) is constrained by the MIDI N-band spectrum and by data from the Spitzer Space Telescope. Our efforts to model 
the structure and SED of IRS 9A range from simple geometrical models of the brightness distribution to one- and two-dimensional 
radiative transfer computations. 

Results. The target is resolved by T-ReCS, with an equivalent (elliptical) Gaussian width of 330 mas by 280 mas (2300 AU by 2000 
AU). Despite this fact, a warm compact unresolved component was detected by MIDI which is possibly associated with the inner 
regions of a flattened dust distribution. Based on our interferometric data, no sign of multiplicity was found on scales between about 
200 AU and 700 AU projected separation. A geometric model consisting of a warm (1000 K) ring (400 AU diameter) and a cool (140 
K) large envelope provides a good fit to the data. No single model fitting all visibility and photometric data could be found, with disk 
models performing better than spherical models. 

Conclusions. While the data are clearly inconsistent with a spherical dust distribution they are insufficient to prove the existence of a 
. disk but rather hint at a more complex dust distribution. 



Key words, techniques: interferometric - stars: circumstellar matter - stars: early-type - stars: formation - stars: pre-main sequence - 
stars: individual: NGC 3603 IRS 9A 

ITl ' 1. Introduction stars we refer to the review articles by McKee & Ostrikeil d2007l) 

" . . . . ., , , and lZinnecker & Yorkd (1200-71). 

H The formation and evolution or massive stars is so rapid that by T , _ , . , , . . . , 

rn 1 ., ^ ( l , i i j .1 • .. • In the case or formation and early evolution of low-mass 

. . , the time the star has reached the zero-age main sequence, it is .„„ .. . . ,. ; TT . . . „ . , 

..„ , , , . , „ , , . , (TTaun) stars and intermediate-mass (Herbig Ae/Be) stars, the 

still enshrouded in large amounts of gas and dust of the molec- , . , • , , „ ■ ■ 

, , , „ , . , . , rrT , • u longer timescales involved allow us to study the mass accre- 

ular cloud from which it was born. Thus, high extinction is usu- . ° , . . , , . . , ; , . , 

„ . . j . V .. , ..r. .. tion process, the chemical evolution of the dust in the regions 

ally present at this stage and prevents the optical identification , , ,, , ■• ■ 

' close to the star, and the subsequent dispersion of the circum- 



stellar material in great detail. A key ingredient of current mod- 



and classification of the star, as well as a more detailed study of 

the inner regions of the circumstellar dust. Furthermore, high- , . 

. , , „ , j. els is the formation of a disk which plays an important role 

mass stars are rare obiects and hence generally much more dis- . . , . r , J . . ,,. r . , 

. . , , . c " , , in removing angular momentum from the inf ailing material, 

tant than low-mass star forming regions, requiring high angular , ■ ■ •, • ^ , j- 

... c . , r , . .? , . £ .i 5 ^ i . Thus, one would expect a similar scenario for the formation 

resolution for study. For a detailed account of the theoretical and . , . ' r . , , , . , . 

, .. , . . j . „ .■ c /u - ■ \ of high-mass stars, even more so since the problem of their 
observational issues connected to the formation of (high-mass) b ... . ■ ■ • , , 
extreme radiation pressure acting upon the spheric al dust en- 

* Based in part on observations collected at the European Southern velopes would limit their fi nal mass (e.g., |Wolfire & Cassinellj 

Observatory, Chile (Prop. No. 074.C-0062) and the Gemini South ll987t lKrumholz et al.l2009l) . But it appears that one cannot sim- 

Observatory, Chile. ply scale the low-mass sta r formation models to higher masses 

** Correspondence: chummel@eso.org dZinnecker & Yorkell2007l) . since, amongst other reasons, most, 
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if not all, high-ma ss stars are born in cluster environments (e.g., 
Ide Wit eta l. 2005) wh ere accretion might occur in a c ompetitive 
manner dBonnell et al.|[T997l 1200 lb IcTark et al.ll2008l) . 

The target of our study, [FPA77] NGC 3603 IRS 9A, has 
been identified as one of the most luminous mid-infr ared sources 
locate d within the giant H II region NGC 3603 dNiirnbergerl 
l200l . This star forming region is dominated by one of the dens- 
est and most massiv e clusters of O and B type main sequence 
stars in the Galaxy (Moffat et al. 1994). The stellar winds and 
ionising radiation of these stars have blown away mos t of the gas 
and d ust of the molecular cloud core around IRS 9A dNiirnbergerl 
2008), which lies at a projected distance of only 2.5 pc from 
the central cluster. In spite of a distance towards NGC 3603 of 
about (7 + l)kpc, the vi sual foreground extinction towards the 
cluste r is only 4-5 mag dSherfl965tlMoffadll983blMelnick et all 
1 19891) . Using infrared colour- magnitude diagram s and pre-main 
sequence evolutionary tracks. iNurnbergen d2003l) was able to de- 
duce a mass of s» 40 Mq for IRS 9A. Its luminosity was esti- 
mated to about 2.3 x 10 5 Lq and its age to ss 10 4 yr. The spectral 
index a^-HVm = 1 -37 is similar to those of (low-mass) IR class 
I objects, indicating that IRS 9A is still embedded in significant 
amounts of gravitationally bound material. 

These properties all argue for IRS 9A's classification as a 
high-mass young stellar object which, due to its advantageous 
location, can be studied more easily. Observations at very high- 
angular resolution in the mid-infrared hold the potential to yield 
valuable constraints for the theoretical understanding of high- 
mass star formation, including multiplicity, but only very few 
other MYSOs have been studied via (mid-) infrared inte rferom- 
etr y so far. Among them are W33A |de Wit et al.ll2007l) . M8E- 
IR dLinz et al.ll2009l) . and MWC 297 dAcke et al.ll2008h . We wiU 
compare the properties of IRS 9A to these objects later on in this 
work. 

The remaining sections of this paper are organised as fol- 
lows. In Sect. |2] we present the main results of our observations, 
shortly discussing the data reduction processes and immediate 
implications of the individual results. In Sect. [3] we model the 
circumstellar structure of IRS 9A using simple geometrical mod- 
els and dedicated radiative transfer models. Finally, we discuss 
the results of these efforts in Sect.|4]and present our conclusions 
in Sect.0 

2. Observations and data reduction 

2.1. MIDI 

Even though the MID-infrared Interferometric instrument 
(MIDI) on the VLTI is primarily an interferometric instrument, 
photometry and imaging are by-products of each observation. 
Thus, a total flux spectrum can be obtained as well as images 
which are taken during the target acquisition procedure. Due to 
the use of an adaptive optics system (MACAO) at the Coude 
foci of each UT of the VLT, these images can be analysed for 
medium scale structures of the target. In the following, we will 
describe each of these three data products. 

2.1.1. Interferometry 

The int erferometric observations using MIDI dLeinert et al.l 
2003b a) were carried out in February and March of 2005. The 
most important details of these observations are summarised in 
TableQ] The resulting (m, v)-coverage is shown in Fig.[TJtogefher 
with the coverage resulting from the T-ReCS observations. The 
synthesized interferometric beam size is about 28 mas by 16 



Table 1. Log of observations of IRS 9A with MIDI, all per- 
formed in 2005. The spatial resolution A/B p is calculated for 
8 \xm. BLPA = baseline position angle projected on sky (east 
of north), AM = airmass. The letters (a)-(d) refer to Fig. [2] 



Date 


Time 


Target 


Bp 


A/B„ 


BLPA 


AM 




(UT) 




[m] 


[mas] 


[°] 




Feb 27 


08h01m 


HD 107446 


38.0 


43.4 


60.6 


1.3 




08h 35m 


IRS 9A 


31.8 


52.0 


79.6 


1.5 




09h 16m 


HD 107446 


33.6 


49.1 


74.0 


1.4 


Feb 28 


06h 24m 


HD 107446 


41.8 


39.5 


44.3 


1.2 


(a) 


06h 48m 


IRS 9A 


38.0 


43.5 


60.6 


1.3 




07h 09m 


HD 107446 


40.2 


41.0 


52.2 


1.2 


(b) 


07h 43m 


IRS 9A 


34.9 


47.3 


70.6 


1.4 




08h 04m 


HD 107446 


37.6 


43.8 


61.8 


1.3 


(c) 


08h 33m 


IRS 9A 


31.6 


52.2 


80.0 


1.5 




08h 54m 


HD 107446 


34.8 


47.5 


70.7 


1.4 


(d) 


09h 21m 


IRS 9A 


28.0 


58.9 


89.9 


1.7 




09h 48m 


HD 107446 


31.0 


53.2 


80.8 


1.5 


Mar 3 


04h 21m 


HD 107446 


56.2 


29.5 


81.4 


1.3 




05h 27m 


HD 107446 


59.4 


27.8 


95.8 


1.3 




07h 10m 


IRS 9A 


62.5 


26.5 


131.9 


1.3 




07h 49m 


HD 107446 


62.5 


26.5 


126.0 


1.3 




08h 54m 


IRS 9A 


61.9 


26.7 


155.5 


1.6 




09h 16m 


HD 107446 


62.2 


26.5 


145.2 


1.5 




09h 33m 


HD 107446 


62.0 


26.6 


149.1 


1.5 




09h51m 


HD 107446 


61.9 


26.7 


153.3 


1.6 



mas, PA of 47 degrees. The observations of NGC 3603 IRS 9A 
were interleaved with measurements of HD 107446, our calibra- 
tor, for which we adopted a diameter of 4. 7 mas based on the 
(V-K) color index (Moz urkewich et al]l2003l) . A star with this 
diameter is nearly unresolved on a 40 m baseline and at a wave- 
length of 8 \im (V 2 = 0.97). HD 107446's angular distance to 
IRS 9 A is about 8°; it has a spectral type of K3.5III and an IRAS 
flux of 32.4 Jy in the 12 \±m filter. MIDI was used with the high- 
sensitivity beam combiner configuration (just two interferomet- 
ric outputs) and the prism for spectral dispersion, providing a 
spectral resolution of R « 35. 

The data were reduced using the MIA+EW S (version 1.5 ) 
software packag43 which is described in detail in Ratzka (120051) . 
In order to calculate the visibility, the correlated flux is nor- 
malised with the total flux. The latter is obtained directly af- 
ter the interferometric measurements, using the individual tele- 
scopes and chopping between the target position and a sky posi- 
tion (for IRS 9A: position angle = 1 17°, throw = 17"). The cal- 
ibration of the target visibilities was achieved by division with 
the calibrator visibilities, interpolated for each epoch of the tar- 
get observations. We estimated the calibration uncertainty from 
the scatter of several measurements of the calibrator on Feb. 28 
to range from about 10% at 13 |^m to about 20% at 8 \xm. The 
resulting calibrated visibilities of IRS 9A are shown in Fig. [2] 
(UT2-UT3 baseline) and in Fig.[3](UT3-UT4 baseline). The sin- 
gle visibility spectrum of Feb. 27 is not shown as it is redundant, 
and the second observation of Mar. 3 was taken at high airmass 
and did not yield data of good quality. The data are made avail- 
able at OLBIN0 

The visibilities are very similar for all our baselines. They 
drop to virtually zero for wavelengths larger than about 9 ii.m 
where IRS 9A appears fully resolved by MIDI. Yet we de- 
tect a compact component (roughly 60 mas in diameter) which 



1 Available at: http : //www . strw.leidenuniv.nl/-nevec/MIDI/index . 

2 http : //olbin . jpl . nasa . gov/data/ 
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Fig. 1. (m, v)-plane coverage of the interferometric observations 
with MIDI at the VLTI and the aperture masking observations 
with T-ReCS at Gemini South. 
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Fig. 2. Squared visibilities measured with MIDI on the UT2- 
UT3 baseline on Feb. 28, 2005, with letters (a)-(d) correspond- 
ing to the four observations listed in Table Q] The solid lines 
are best fits of the geometrical model described in Section I3.ll 
and the dashed lines are from the physical model described in 
SectionEHl 



emerges below this wavelength, causing a steep rise of the visi- 
bilities. 



2.1.2. Spectro-photometry 

Since HD 107446 is not a spectro-photometric calibrator, we 
adopted for its spectrum the template of HD 16 3588 (K2III) 
from t he list of spectro-photometric standards by ICohen et al.l 
( 1999!), and used the ratio of the IRAS 12 \im fluxes to deter- 
mine the relative scale. The resulting flux-calibrated spectrum of 
IRS 9A as measured with MIDI is shown in Fig. [4] We averaged 
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Fig. 3. Squared visibilities measured with MIDI on the UT3- 
UT4 baseline (first observation on Mar 3, 2005). The solid line 
is a best fit obtained from the geometrical model described in 
Section [3~Tl and the dashed line is from the physical model de- 
scribed in Section 1331 
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Fig. 4. Average spectrum of IRS 9A, as obtained from the MIDI 
data taken on Feb 28, 2005. The small spike at 12.8 [im coincides 
with the [Nell] line. 



only spectra from Feb. 28, 2005 as they represent the largest ho- 
mogeneous data set. The error bars are based on the standard 
deviation of the spectra. 

There are several interesting aspects concerning this spec- 
trum. First, the flux of IRS 9A rises steeply with increasing 
wavelength. Second, there is no sign of a silicate feature at about 
9.7 |xm, neither in absorption nor in emission. This is rather 
uncommon for dust-enshrouded objects like YSOs. Third, de- 
spite the low spectral resolution of MIDI, there is a hint of the 
[Nell] emission line at 12.8 [im, prominently seen by Spitzer 
(see Sect. I2.31>. 



2.1.3. Imaging 

The chopped acquisition images of IRS 9A were taken with the 
N8.7 filter (A c = 8.64 \xm, /Ifwhm = 1-54 (j.m) and have a spatial 
resolution of about 0.3". Using the calibrator HD 107446 as a 
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Table 2. Results of the Gaussian fits to the deconvolved MIDI 
acquisition images. 
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Fig. 5. Deconvolved MIDI acquisition image of IRS 9A in the 
N8.7 filter after 25 iterations. The contour lines are equally 
spaced for the fourth root of the intensity and lie at 3.7, 6.3, 10, 
15, 22, 32, 44, 59 and 77 percent of the maximum intensity. 



PSF reference, the images can be deconvolved to further enhance 
the res olution, using the procedure described bv lChesneau et al.l 
(2005). We used the Lucy-Richardson algorithm (Luc vTl974l) for 
this purpose, taken from the IDL Astronomy Library^ However, 
we detected a residual jitter in our images caused by a delay in 
the acquisition by the adaptive optics loop each time the chop 
position changed. Simply re-centering the frames did not lead to 
useful results due to distorted profiles while MACAO was clos- 
ing the loop. Therefore, we only stacked the chopped frames 
which were well centered (about half of the frames). 

We stopped the algorithm after 25 iterations, a compromise 
between image convergence and fitting noise. A representative 
example of our deconvolved acquisition images is shown in 
Fig. [5] We performed two-dimensional Gaussian fits to the de- 
convolved images in order to determine the FWHM sizes and 
the orientation of the major axis. The results are summarised in 
Table [2] which lists the mean single night values, and also the 
mean value of all our observations, using the number of obser- 
vations per night as the associated weight (N t, s ). We did not use 
the images of the beam B of MIDI (UT3 on Feb. 27 and Feb. 
28, UT4 on Mar. 3) due to their slightly inferior image quality 
and the associated systematically larger FWHM values. Overall, 
IRS9A has a Gaussian FWHM of about 389 mas by 354 mas, 
and a position angle of the major axis of about 75 degrees (east 
of north). Despite the different values of the parallactic angle of 
the observations, i.e. the direction to north on the detector, rang- 
ing from 130 to 173 degrees, the derived PA of the fitted ellipse 
was unchanged. Representing the overall shape of the source, 
the fitting of the PA was apparently dominated by the SW ridge 
in the source structure. 



2.2. T-ReCS 

We also observed NGC 3603 IR S 9A with the T-ReCS (The rmal- 
Region Camera Spectrograph, Ide Buizer & Fisherl (l2005h ~) in- 
strument at the Gemini South Observatory on January 2, 2005. 
A 7-hole non-redundant aperture mask was used to reduce the 
weight of the short spacings realised within the 8 m aperture of 
the telescope, and therefore to yield high quality images at the 



Date 


Beam 


N obs 


Major 


Minor 


PA 








[mas] 


[mas] 


n 


Feb. 27, 2005 


UT2 


1 


402 


339 


82 


Feb. 28, 2005 


UT2 


4 


392 


366 


71 


Mar. 03, 2005 


UT3 


2 


375 


339 


81 


Mean values 




7 


389 


354 


75 



1.0 
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T3 



cr 
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Projected baseline length [Ml] 
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Fig. 6. Squared visibility amplitudes from T-ReCS plotted ver- 
sus projected baseline length. Error bars have been omitted for 
clarity; they are shown in Fig. [7] The vertical scatter of the visi- 
bility amplitudes arises from non-axisymmetric structure, which 
causes some azimuthal dependence of the amplitudes. The solid 
line is a best fit obtained from the geometrical model described 
in Section I3.ll and the diamonds correspond to the physical 
model described in Sectionl3~3l 



diffraction limi t (for a description of the aperture masking tech- 
nique see, e.g.. lMonnierll 999: Tut hill et al.l2 000). The data were 
taken with the Si-5 filter (A c + aA = 11.66 |am+0.57 |am). The re- 
sulting (m, v)-coverage is shown in Fig.[T] the synthesized beam 
is almost circular with a diameter of 264 mas. 

The interferograms, which result from the superposition of 
all interference patterns between any combination of apertures in 
the mask, are Fourier transformed to measure the visibility am- 
plitudes and phases. This is done for both target and calibrator, 
and the calibration procedure is similar to the one used in long 
baseline interferometry. The resulting squared visibility ampli- 
tudes are shown in Figs.[6]and[7] The data are made available at 
OLBIN0 

In order to derive an image from the set of amplitude s and 
phase s, we employed the VLBMEM imaging software dSivial 
1987), which uses a Maximum Entropy algorithm. The resulting 
image is shown in Fig. [8] It agrees well with the overall appear- 
ance of our MIDI acquisition images, again showing the emis- 
sion of IRS 9 A to be quite extended and slightly asymmetric. 



3 http://idlastro.gsfc.nasa.gov/ 



4 http : //olbin . jpl . nasa . gov/data/ 
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Fig. 7. Squared visibility amplitudes from T-ReCS (+- symbols) 
plotted versus baseline position angle (0-180) The vertical off- 
sets of the visibility amplitudes is due to the different lengths of 
baselines with the same position angle. The solid lines (with x- 
symbols) are the best fits obtained from the geometrical model 
described in Section[3~Tl 
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Fig. 8. T-ReCS image of IRS9A at 11.7 nm. North is up, and 
east is to the left. Reconstruction from the aperture masking data 
using the Maximum Entropy Method (MEM). Contour levels are 
the same as in Fig. [5] 



2.3. Spitzer 

In Fig. |9] we show the SEP of IRS 9 A measured by 
iLebouteiller et al.l(l2007ll2008h with the Spitzer Space Telescope 
in comparison with the MIDI spectrum. Since Spitzer's spa- 
tial resolution is about a factor of 10 worse compared to the 
(MACAO-assisted) VLT, the IRS spectrograph aboard Spitzer 
collects considerably more flux, about twice that of MIDI in the 
middle of the N band. There is no hint of the [S IV] line in our 
MIDI data, although its line flux is m ore than two times highe r 
than the one from the [Ne II] line (see lLebouteiller et alJ (|2008)) 
which we did barely detect with MIDI (see Fig.|4]i. It is possible 
that the region emitting the [S IV] line is more extended than the 
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Fig. 9. The spectrum of IRS 9A from Spitzer (solid lines) and 
MIDI (dotted line). The numerous forbidden emission lines are 
labelled together with their central wavelengths. There are also 
several PAH features visible, with the most prominent ones ly- 
ing at 6.2, 7.6 and 1 1.2 [im. The small gap is due to saturation of 
Spitzer's SL module for wavelengths larger than 9.7 |am. The 
data from the LH module (18.7-37. 2 |i.m) suffered from th e 
same effect and are not shown (see also Lebouteil ler et al.l2 008). 



MIDI field of view. We will use the Spitzer SED as an additional 
constraint for our models in Sect. [3] 



3. Model fitting and analysis 

In this section, we first try to model the circumstellar structure 
of NGC 3603 IRS 9A in the context of current star formation 
scenarios by using such elements as rings, disks, and envelopes. 
We then describe our efforts of testing physical models which 
require radiative transfer computations to produce images of the 
structure for comparison with the available visibility and photo- 
metric data. However, given that the environment of young and 
very luminous stars is complex, our models can only represent a 
first step on the way of developing an adequate physical model 
of IRS 9A. Also the paramerters of these models may not be well 
constrained given their large number. 



3.1. Simple geometrical models 

It is useful to identify structural components of the source and 
derive the parameters of such simple geometrical models like 
sizes and orientations. Beginning with the visibilities measured 
with T-ReCS in a single filter (hence no wavelength dependence 
of the source structure need to be considered), we first fit circu- 
lar Gaussian and uniform disk models with diameters of 3 1 1 mas 
(FWHM) and 495 mas, respectively. The reduced^ 2 (henceforth 
denoted with Xr) of the Gaussian model can be further reduced 
to 2.3 by allowing an ellipsoid, with a major axis of 329 mas, 
PA of 100° east of north, and a minor axis of 282 mas (axial ra- 
tio 0.86, compare with values from Table [2}. This component 
would be completely resolved by MIDI even on the shortest 
baselines available, and therefore a second component of small 
angular scale needs to be added to the geometrical model in or- 
der to explain the MIDI detection. After trying unsuccessfully 
to fit the data with an additional uniform disk, we adopt a ring- 
shaped component instead because its Fourier transform exhibits 
a larger secondary maximum (beyond the null) than a uniform 
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Table 3. The parameters of the geometrical model. 



Parameter 


Unit 


Value 


Ring major axis 


[mas] 


57 


Major axis PA 


[°] 


105" 


Ratio minor/major axis 




0.32 


Fractional ring width 




0.31 


Ring effective temperature 


[K] 


1000 


Envelope major axis FWHM 


[mas] 


350 


Major axis PA 


[°] 


105" 


Ratio minor/major axis 




0.86 


Envelope eff. temperature 


[K] 


140 



Forced to be the same value. 



disk. The ring, which could be associated with the inner hot re- 
gions of a disk with a central hole, is initially assigned a black- 
body temperature of 1500 K (equal to the typical dust sublima- 
tion temperature) and a fractional width of 1/4 of its radius. We 
assume here, before a more realistic physical model is applied 
later on, that the two components are co-located and do not ob- 
scure each other. In order to have the flux of the Gaussian com- 
ponent decrease towards the blue end of the N-band so that the 
correlated flux of the ring begins to dominate here the visibil- 
ity measured by MIDI (see Fig. ITOb . we assign initially a low 
blackbody temperature of 250 K to the exte nded Gaussian com - 
ponent based on a spectral decomposition b v lNiirnbergerl (12003 ) . 
This model ("geometrical model") indeed results in a good fit 
(Xf = 1.6) of all visibility data (see Figs. [2] and [3] as well as 
Figs.|6]and|7]i, with the best-fit parameter values given in Table[3] 
The inclination of the polar axis of the ring against the line of 
sight is 71 degrees. The fit to the SED is shown in Fig.lTOl 

There is no clear evidence for a companion of IRS 9A, and 
we can use the geometrical model to estimate at what level we 
can exclude a companion. We adopt the companion to be a ring 
like the primary so that the binary model asymptotically ap- 
proaches the geometrical model for small separations when the 
individual component fluxes are adjusted (given a magnitude dif- 
ference) to yield a constant total flux. For the limiting xf> we 
adopted 1.7 which results in recognizable deviations from the 
data. Almost all models below this limit have a magnitude dif- 
ference of 3 or larger, and therefore we can exclude the possi- 
bility that a companion of more than about 0.5 Jy exists between 
30 mas (the resolution limit) and 500 mas (our search area) away 
from the primary. Due to the fact that uncorrelated flux domi- 
nates at the longer wavelengths where none of the models could 
therefore produce any significant signatures, we cannot place 
more stringent limits. 

3.2. Spherical dust shell models 

We now interpret the geometrical model we have studied as 
spherical dust shells with the innermost one corresponding to 
the ring component and with the star at their center. For the 
calculation of the radiative transfer we use the DUSTY code 
(llvezic&Elitzurlll997l 119991: llvezicet al.lfl999h which adopts 
spherical symmetry. It produces detailed spectra and images 
at specified wavelengths, and the latter can be used to calcu- 
late the visibilities of the models that correspond to the ob- 
servational setup of T-ReCS or MIDI . These type s of m odels 
have been u s ed qu ite successfully by Ide Wit et alj d2007l) and 
Ide Wit et al.l d2009h for their MIDI observations of W33A and 
resolved 24.5 \xm emission of MYSOs, respectively. 

For the central star we adopted an effective temperature of 
22 000 K, and for the shells a density, p, proportional to r~ 2 , 




1 10 100 

Wavelength X [urn] 



Fig. 10. The SED of the geometrical model (solid line) with the 
observations from Spitzer (green and blue solid lines), MIDI (red 
sol id line), and NIR/MIR photometry (plus symbols) published 
by iNiirnberperl (120031) . The two components of the model have 
the same flux at 8 microns. The dashed and dotted lines corre- 
spond to the physical model described in Section [3~3l with aper- 
ture radii of 3" (Spitzer) and 0.3" (MIDI), respectively. 

where r is radius of the shell. For the dust composition we used 
standard astronomical silicates. The most important free param- 
eters to adjust were the temperature, T s , of the inner edge of 
the shells, and its opacity, r. In general, a lower T s increased 
the size of the shells and thus lowered the visibility amplitudes, 
while a larger opacity had a similar effect yet introduced a strong 
silicate absorption feature in the spectrum. Using the T-ReCS 
data and the SED as constraints, we settled on a model with 
r s = 500 K and Tip„ m - 1.5, corresponding to Ay - 29 mag 
(Rie ke & Lebofskvlil985l ). The fit to the T-ReCS visibilities was 
indistinguishable from the (circular) Gaussian model, while the 
fit to the SED is shown in Fig. Q~T] DUSTY returns the radius 
of the innermost shell, r m (where T = T s ), which scales propor- 
tional with the square root of the luminosity, L 1 ^ 2 . Adopting a 
luminosity of IRS9A of 230 000 Lq, we derive r s = 170 mas. 
The (dust) mass of the shell can be derived from the opacity r, 
the absorption coefficient k, and the ratio of the outer edge of the 
shell to r m , Y, as follows, 

r =£ m Kpds =p K— (r out -%), (1) 

,n f out 



M = pAnr 2 dr - 4n-rf n Y, (2) 

giving a dust mass of about IMq. 

We thus find that it is possible to devise a physical model 
for the T-ReCS data set, but when we apply that model to the 
observations from MIDI, the agreement is very poor because the 
emission is completely resolved spatially. This also applies vice- 
versa, that is we are able to find models with DUSTY which 
can reproduce the visibilities measured by MIDI (albeit showing 
larger deviations), but these models are in turn a very poor fit 
to the data from T-ReCS. In order to find suitable models that 
are able to account for all our observational data, we conducted 
a grid search varying all parameters except for T s , Y and the 
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E 

* 




10 

Wavelength A [wm] 



100 



Table 4. The parameters of the Robitaille disk-envelope model 
3012790. 



Parameter 


Unit 


Value 


Stellar mass 


[Mq] 


25 


Stellar radius 


[*©] 


6.5 


Effective temperature 


[K] 


38000 


Luminosity 


[L©] 


92000 


Inner disk/envelope radius 


[AU] 


25 


Outer disk radius 


[AU] 


94 


Outer envelope radius 


[AU] 


100000 


Disk dust mass 


[Mq] 


0.005 


Envelope dust mass 


[Mq] 


0.9 


Inclination 


n 


85 


Disk flaring power, /? 




1.2 


Disk scale height 


[AU] 


9 


Cavity cone angle 


[°] 


29 



Fig. 11. Comparison of the SED of our DUSTY model with the 
observations from Spitzer, MIDI, and the NIR/MIR photometry 
(plus symbols) from Fig. [10] The dashed line corresponds to the 
attenuated input radiation, the dash-dotted line to the dust radia- 
tion, the dash-double-dotted line to the scattered emission. 



temperature of the central star. However, none of the resulting 
models could accomplish this task. 

It thus seems obvious that spherical dust distributions are 
not compatible with the circumstellar structure of IRS 9A. 
Therefore, in the next section, we examine the effects of non- 
spherical dust distributions on the visibilities and the SED. 

3.3. Disk models 

Within the framework of models for YSOs, the most obvious 
approach for a non-spherical density distribution is given by the 
flattened structure of a circumstellar disk embedded in an enve- 
lope. If seen under low (near face-on) or intermediate inclina- 
tions of the disk axis to the line of sight, a disk opens a more 
or less direct optical path to the central object and the hot in- 
ner disk regions (the hot ring in our geometrical model), while 
at the same time maintaining an extended dust distribution (the 
cool Gaussian component of the geometrical model). Models 
involving disks have already been used widely for interpreting 
SEDs and interferometric observations of T Tauri and Herbig 
Ae/Be stars. For example, the visibilit ies predicted by m odels 
computed for Herbig Ae/Be stars by iLeinert et alj (|2004) look 
similar to our MIDI visibilities. In this section we report on our 

first attempt to fit such models to our data. 

We adopt the disk-en velope model of Whitney et al.l (|2003b) 



which has been used by Robitaille et al. (2006) to compute a 
large grid of model SEP, ma de available through a web servefl 
Following iLinz et all ([2009), we select only those models fit- 
ting the Spitzer-SED of IRS9A and compute their images. To 
compute the model visibilities, we apodized the images at scales 
larger than the PSF of an 8 m telescope in the case of MIDI, 
while due to the comparitively small subapertures in T-ReCS, 
this apodization was not necessary. Thus we find that the model 
3012790 inclined at 85° to the line of sight (nearly edge on) and 
with the disk roughly oriented East- West provides the best fit of 
the ten selected models. The parameters of this model are listed 
in Table|4] The disk-envelope model 3012790 includes a signifi- 



http : //caravan . astro .wise . edu/protostars/ 



cant accretion rate, for which however no independent constraint 
exists so far. The fits of this model to the visibilities are shown 
in Figs.[2][3j and Fig. [6] while the SED is shown in Fig.fTUl 

4. Discussion 

Considering the MIDI observations, the most important feature 
to reproduce by any model is the steep rise of the visibility to- 
wards 8 microns, indicating the emergence of warm inner dusty 
regions. The disk-envelope model is able to reproduce this fea- 
ture about as well as the geometrical model, but the envelope 
is too large considering the low predicted visibilities for the T- 
ReCS observations. 

The T-ReCS image of IRS 9A (Fig. [8) rese mbles the syn- 
thetic i mages of class I sources computed by IWhitnev et al.l 
(2003a) (their Fig. 11a) and indicates the presence of a bipo- 
lar cavity facing the observer and pointing towards the South to 
South- West. Ionizing stellar radiation escaping through the cav- 
ity has created an HII region as indicated by the detection of the 
[Nell] line by Spitzer (see Fig. |9j. This morphology of the cir- 
cumstellar dust (and gas) distribution is consistent wit h the inter- 
pretati on of NACO images of NGC 3603 IRS 9A bv lNurnbergeri 
(2008) which show an envelope with a major-to-minor axis ratio 
of about 2:1 at a position angle of 82°. (The NACO observa- 
tions also confirm the reality of the faint feature in our T-ReCS 
image about 0.6 arcseconds to the north of IRS 9A.) Due to the 
limited dynamic range of our images, we cannot confirm the ab- 
lation of dust and gas towards the south-east, i.e. away from the 
NGC 3603 cluster center, as seen in the NACO images. 

A lower mass estimate of 0. 1 Mp) for the hot du st and gas 
bound to IRS9A was derived by iNiirnbergerl (120031) based on 
the assumption that the MIR emission is optically thin. For 
the envelope mass, an estimate of 1 Mq is derived by the 
same author based on an estimate of the intrinsic extinction of 
Ay = 10 -15, much lower t han our estimate based on DUSTY. 
However, Niir nbergerl d2003l) remarked that a radio flux measure- 
ment of NGC 3603 IRS 9A yielded unexpectedly high values, 
which would argue in favour of a higher envelope mass estimate. 
A final answer cannot be given until a sub-mm flux measurement 
can be obtained. 



5. Conclusions 

We gave a detailed account of our efforts to model the observed 
spatial and spectral data of the MYSO candidate NGC 3603 
IRS 9A. As a first result, we find that the observed visibilities are 



8 



S. Vehoff et al.: Mid-infrared interferometry of NGC 3603 IRS 9A 



not compatible with a binary nature of IRS 9A. However, given 
the limits in sensitivity (» 0.5 Jy; the correlated flux on longer 
baselines would have been too small to detect) and in angular 
resolution (down to 30 mas) of our interferometric data, we can- 
not exclude the presence of a companion with either significantly 
lower mid-infrared flux (less than 10%) and /or an extremely 
tight binary with a separation of less than about 100 AU. 

We have shown that a spherical dust distribution is incon- 
sistent with the data and that therefore a flattened distribution 
is more likely. A geometrical model consisting of a warm ring 
and a cool extended component can reproduce all observations 
fairly well, but requires the ad-hoc assumption of a direct line 
of sight into the inner regions of the warm dust. The ring has 
a temperature 500 K less than the dust sublimation temperature 
(1500 K), probably a consequence of the fact that it represents 
not only the very inner edge of a disk, but rather the inner regions 
of it. Physical disk-envelope models provide a self-consistent 
connection between the hot inner and the cooler outer regions 
and allow the former to be seen more directly. However, no sin- 
gle model was able to fit all observations simultaneously. Since 
it could be expected from the early evolutionary stage of class 
I sources compared to the well developed disks of the class II 
sources (adopting the classification scheme more strictly devel- 
oped for lower mass objects), that the dust structures are more 
complex, e.g. clumpy or showing disk distortions, our data are 
not sufficient to prove the existence of a disk. 

Comparing NGC 3603 IRS 9A to other MYS Os for which 



MIDI interferometry was published, e.g. W33A (de Wit et al 

2007) , M8E-IR dLinz et all 12009). and MWC 297 lAcke et al 

2008) , we find that all of them exhibit rather different properties. 
These depend mostly on the amount of circumstellar material 
present, as reflected in the estimated visual extinctions which 
range from Ay = 22 to Ay = 100. The source with the high- 
est extinction, W33A, displays a deep silicate absorption fea- 
ture. The measured visibilities were similar to the ones we mea- 
sured between 8 and 9 i^m on the same baseline, translating into 
a Gaussian FWHM of 120 AU for W33A (distance of 3.8 kpc). 
A significant difference here to NGC 3603 IRS 9A is that the lat- 
ter is completely resolved by MIDI at longer wavelengths due 
to the envelope detected with T-ReCS. Thus, W33A seems a lot 
more compact, and was therefore successfully modelled with a 
single DUSTY envelope. The model required a rather low effec- 
tive stellar temperature for a 10 Mq star (only 10 000 K), possi- 
bly indicative of a star swollen due to accretion. The same effect 
was observed in M8E-IR, which was modelled with an envelope 
and disk surrounding a rather cool (5 000 K) star with a mass 
of 13.5 Mq. M8E-IR has a much weaker silicate absorption fea- 
ture, and appears to be rather compact (45 - 100 AU). The ex- 
istence of a disk in the envelope could not be ascertainted. Yet 
again, the difference to NGC 3603 IRS 9A lies in the fact that the 
visibility of M8E-IR is only weakly dependent on wavelength, 
indicating that a flattened structure is not necessary to explain 
them. Finally, the observations of MWC 297 could not be mod- 
eled with a physical disk model at all. This star is a young 10 Mq 
Be star with a rather low extinction. Only indirect arguments 
favoured the disk geometry for the circumstellar material. 
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